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ABSTRACT 

We consider a general class of boundary terms of the open XYZ spin- 1/2 chain 
compatible with integrability. We have obtained the general elliptic solution of 
i^-matrix obeying the boundary Yang-Baxter equation using the i?-matrix of the 
eight vertex model and derived the associated integrable spin-chain Hamiltonian. 



1+1 dimensional integrable models with boundaries find interesting applica- 
tions in particle physics as well as condensed matter systems. In view of this, 
attempts have been made recently at the integrable extension of conformal field 
theories, ' both massive and massless integrable quantum field theories and 
solvable lattice models' 7 " 1,31 to those with boundary terms. In the case of lattice 
models, relying on the earlier work by Cherednik,' 81 Sklyanin has given 191 a general 
framework which enables us to treat this problem on an algebraic footing. Es- 
pecially, the general solution of integrable boundary terms has been found in the 

m 21 

XXZ and the XXX Heisenberg spin chain system based on his framework. In 
this letter, we will work out the general solution of boundary interactions in the 
case of XYZ Heisenberg spin chain system. 

The Hamiltonian of the XYZ spin- 1/2 chain is given by the transfer matrix 
of the eight vertex model.' 141 The eight vertex model is defined in terms of the 
Boltzmann weights given by the elliptic solution R(u) of the Yang-Baxter (Y-B) 
equation 



Ru(u - u')R 13 {u)R 2S {u') =R 23 (u')R 13 (u)R 12 (u - u'). 



(1) 



Here we regard R{u) as linear operators acting on the tensor product of vector 
spaces V <g> V with V = Cv + © Cv _ and R\ 2 = R <S> 1, -R23 = 1 <8> -R, etc. as 
those acting on V\ <g> V 2 <g> V3, where Vi = V, i = 1, 2, 3. Setting R^v^ (g> v £ > 2 = 
^2 £l£2 v ei <H> v £2 R{u) £ e \el and arranging the elements of R in the order {e\,e 2 ) = 
(++), (H — ), ( — h), ( ), one can express the eight vertex i?-matrix as follows. 



R{u) 



( sn(w + rj) 






k sn i] sn u sn(w + rj) \ 
sn u sn i] 
sn rj sn u 



\ k sn rj sn u sn(w + rj) 



sn(u + rj) 



(2) 



/ 



where sn u = sn(w; k) is the Jacobi elliptic function of modulus < k < 1. 
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Let Vij be the transposition operator on V% ® Vj, i.e. V(x <S> y) = y <E> x. The 
i?-matrix (2) is known to have the following desirable properties. 

Regularity R(0) = r{rj)V, r(rf) = sn rj, (3) 

P- symmetry V\<iRvi{u)V\2 = Ru(u), (4) 

T - symmetry R\ 2 t2 (u) = Ru(u), (5) 

Unitarity Ri2(u)Ru{—u) = p(u)l, p{u) = sn 2 rj — sn 2 u, (6) 
Crossing unitarity R^u) R^i—u — rj) — p(u)l, 

p{u) = sn 2 -q — sn 2 (u + rj). (7) 



In the case of periodic boundary condition, it is known that the Y-B equation 
(1) implies a commuting family of transfer matrix. Hence the model is integrable. 

We now consider the eight vertex model with boundary interactions. Aiming 
at describing integrable systems with boundaries, Sklyanin' 91 has introduced a pair 
of matrices K+(u) and K-(u). The effects of presence of boundaries at the left 
and right ends are solely described by K + (u) and K-(u), respectively. K±(u) are 
defined as the solutions to the relations 

R 12 (u - u) K- (u)R 12 {u + v!) K- (u) 

=K- (u')R 12 (u + u) A_ (u)R 12 (u - u'), (8) 

1 2 

R 12 (-u + u') K% (u)R l2 (-u - u' - 27]) K% («') 

2 1 

=K% ( U ')R 12 (-u - u' - 2ri) K% (u)R 12 (-u + u'), (9) 

1 2 

where K±= K± ® idy 2 and K±= idy 1 <E> K±. The equations (8) and (9) are called 
boundary Y-B equations and K±(u) boundary iC-matrices. 

The boundary Y-B equations imply a commuting family of transfer matrix.' 91 
The transfer matrix t(u), in this case, is defined using the K± and the monodromy 
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matrix T{u) as 



t(u) =tr K+^T^K^T'^-u) , (10) 
where T(u) is given by 

T(u) =R N0 (u)R N - W (u) ■ ■ -Rw(u). (11) 
The trace in (10) should be taken over Vq. Then, the commuting property of t(u), 

[t(u), t(u')\ =0, (12) 

follows from the properties of R and the boundary Y-B equations (8) and (9). 

The problem is now to solve the equations (8) and (9) and find general solutions 
for K- and K+, using the eight vertex i?-matrix given in (2). It suffices to consider 
the first equation, because of the following fact. Suppose K-(u) is a solution of 
the first equation, then the function 

K+(u) =Ki(-u -rj) (13) 

gives the solution for the second equation. 

We now proceed to solving Eq.(8). Write K-(u) as 

K-{u) = , (14) 

\z(u) v(u) J 

we have found that, out of the sixteen equations in the boundary Y-B equation 
(8), only three are independent: 

S-VV 1 + S + XV f =S + Vx' + S-XX* , (15) 

yz + k s-S+zz' =zy' + k s-S+yy', (16) 



S-s + yx' + k s^S-S-zx' + k s v S-S+(s-vz' + s+xz') 

=S + s-yx' + k s^S+s+zx' + s v (s-vy' + s+xy'), (17) 

where we set x = x(u), x' = x(u'), etc., and 

s v =sn rj, s± = sn(u±u), S± = sn(u ± u + rj). (18) 

In the following, we also use the notations a(u) = v(u)/x(u), f3(u) = z(u)/y(u) 
and j(u) = y(u)/x(u). 

Dividing (15) by xx', one obtains 



sn(ii + u')a(u) + sn(u — u') 

a \ u ) — — 7 K — 7~^ 7 77- (19) 

sn(w — u'Ja(tij + sn(w + w) 



Taking the limit u' — ► u of Q ^J_^^ , one obtains the following differential equa- 
tion. 

da ( u ) = _ l-a(u) 2 
du sn 2u 

After the change of variable t = sn u, the integration of (20) takes the form 

f da 1 f l-k 2 t 4 

J - - 2 J d t(i-t 2 )(i-m 2 y { } 

One can easily get the general solution of Eq.(21), 

v(u) Ccn u dn u — sn u 



x(u) Ccn u dn u -\- sn u ' 
where C is an arbitrary constant. 



(22) 



In a similar way, from (16) one gets 



, _ P(u) + k sn(u + u')sn(u - u') 
P{U> k sn(u + u')sn(u - u')(5(u) + 1 ' 1 ' 



and the differential equation 

<W --ksn2u(l-(3(u) 2 ). (24) 



dw 

This implies the general solution 

z(u) A(l — fc sn 2 u) — 1 — A; sn 2 u 
y(u) A(l — fc sn 2 tt) + 1 + fc sn 2 u' 

with A being another arbitrary constant. 

Dividing (17) by j/y', the third equation (17) can be written 

7(u) (s_a(u) + s+)(l — k s + s_/5(-u / )) 



(25) 



7(m ; ) - S+s- + fc /3(m)s 2 (S-S- - S+s+) ' 



(26) 



Substituting (23) into (26) and replacing ct(u) and f3(u) by the RHS of (22) and 
(25), one can factorize (26) in the form of the ratio of the same functions, one with 
the argument u and the other with u', respectively. We thus find 

y(u) _ A(l - k sn 2 u) + 1 + k sn 2 u 
x(u) Ccn u dn u + sn u 

where \i is the third arbitrary constant. 

From (25) and (27), one can now get the ratio z(u)/x(u). It is then not difficult 
to check that the above solutions satisfy all the remaining equations, if one notes 
the identity 

S ^ Ss + - S + s- = l-k*S + Ss + s- = S « Sn(2M + V) - (28) 
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In summary, we have obtained the general solution of (8) as K-(u) = K(u; £, A, 
with 



/ (j- , \ r, A(l— k sn 2 u)+l+fc sn 2 it 

_L_ / sn(£ + «) /i sn2w an , u 

„ n f n A(l— k sn 2 u) — 1— k sn 2 u \ 
sn 4 V /i ™2u 13^%^ sn(£ - «) 



(29) 



where we set C = sn£/cn£ dn£ and replaced jj, cn£ dn£ by ji. We normalize the 
matrix K-(u) as if-(O) = 1 for later convenience. 191 

In the trigonometric limit k — > 0, where sn u — > sin u, we recover the result in 

/ [12] 

the case of the six vertex model given by de Vega and Gonzalez Ruiz. The rational 
limit is obtained from the trigonometric i^-matrix by scaling u — > rju, £ — > 77^ and 
taking the limit 77 — ^ 0. 

Let us next consider the corresponding XYZ spin chain Hamiltonian. Because 
of the equation (12), one can regard the transfer matrix t(u) as the generating 
function of integrals of motion of the system. Its first logarithmic derivative implies 
the Hamiltonian. 

H =2r(r]) ^ 1 (0)(t / (0) - tri^(0)) 

^ig^^-Wl ^* ), (30) 



71=1 

where the two-site Hamiltonian is given by 

H n ,n+1 =^-jVnn+lRnn+l(fy- ( 31 ) 

By a direct calculation with the K-matrices K-(u) = K-(u;^-,X-,/i-) and 
K+(u) = 



-k The solution obtained in Ref.[12] for the XYZ model are the special cases of (29) associated 
with the special solutions a(u) 2 = 1, (3(u) 2 = 1 of Eqs.(20) and (24). 
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KL(—u—r); — £+, — A+, —(J,+) together with the i?-matrix (2), one gets the following 
result. 



N-l 



h = Yl ((! + r )«+i + (1 " r )*I + i + A« +1 ) 



n=l 



+ sn 77(A_of + 5^(7+ + CL^ + + B + a+ + C+o"^), (32) 



where 



T =/c sn ?7, A = cn rj dn rj 

cn £ ± dn g± /x ± (A ± + 1) /J±(A± - 1) , , 

2sn ^ ± sn i ± sn ^ ± 

In conclusion, we have obtained the general elliptic solution of the boundary 
Y-B equation for the .fT-matrices and derived the Hamiltonian of the associated 
XYZ spin- 1/2 chain with boundary terms. 

An immediate question is to find the ground state energy and the excitation 
spectrum of the XYZ Hamiltonian we have derived. The diagonalization of this 
Hamiltonian can be achieved by means of the generalized algebraic Bethe ansatz by 
Sklyanin' 91 with some modification as in the periodic boundary condition case.' 14 ' 151 
This subject is now under investigation. 

It was shown using the results of Bethe ansatz type analysis that, by tuning 
the X-Y anisotropy coupling (T in Eq.(32)), the XYZ Hamiltonian with periodic 
boundary condition give rise to the quantum sine-Gordon theory in the continuum 
limit! 16 ' In the case of open XYZ spin chain, it is of interest to ask whether one can 
tune the coupling of boundary terms together with the T so that one can derive its 
field theory limit. The resulting theory is expected to be the quantum boundary 
sine-Gordon theory' 3 ' 41 In the limit iV — > oo, we have three boundary terms pro- 
portional to af, o\ and erf associated with three free parameters in K_, whereas 

T171 [31 

the boundary term proposed by Sklyanin and Ghoshal and Zamolodchikov 
has two parameters. It is necessary to explain this difference in the analysis of the 
continuum limit. 
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Furthermore, the higher logarithmic derivatives of commuting transfer matrix 
give infinite number of conserved quantities. In the closed XYZ spin chain case, 

[18] 

it was shown that the conservation laws associated with these quantities yield 
selection rules in the scattering process of the quantum sine-Gordon theory. It is 
an interesting question to ask how the parameters appearing in the i^-matrices 
affect the scattering process of the boundary sine-Gordon theory. 

We will present our result on these problems in future publications. 
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